Hybridization of localized and density modes for the roton spectrum of superfluid 4 He. 
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A new description is presented for the roton part of the energy spectrum of superfluid 4 He. It is based 
on the assumption that in addition to the Feynman density excitations, there exist localized modes 
which describe vortex-core elements. The energy spectrum which results from the hybridization of 
these two kinds of excitations is compared to the experimental data namely the structure factor 
and the scattering cross-section for the single quasiparticle excitation. Another type of excitation 
interpreted as a vortex loop is obtained whose energy agrees both with Raman scattering and critical 
velocity experiments. 
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In this letter we present a new theoretical descrip- 
tion of the roton region of the excitation spectrum of 
superfluid i He. In spite of the success ^tained using 



the wavefunction proposed by Feynman un, a quantita- 
tive agreement with the experimental results has not yet 
been achieved. Subsequent refinements using more com- 
plicated variational wavefunctions were at the expenses 
of a simple interpretation as a pure density, fluctuation^. 
More recently, new experimental results □ obtained in 
high resolution neutron scattering did show puzzling fea- 
tures such as the persistence of the phonon peak up to 
temperatures higher than T\ together with the decrease 
of the sharp peaks associated respectively to the maxon 
and roton excitations. These features tend to indicate 
that the low momentum phonon mode is marginally re- 
lated to superfluidity or even to the Bose- Einstein statis- 
tics, a point which was already noticed by Pines u who 
proposed to interpret this mode as a collective zero sound 
just like in 3 He. In contrast, the decrease at the transi- 
tion of the other modes suggests a connection with quan- 
tum exchange effects. 

In order to interpret these features we assume the ex- 
istence of two kinds of excitations in the system. At 
large scales, there are delocalized density fluctuations 
while at scales of a few interatomic distances there ex- 
ist localized excitations which may be associated with 
microscopic vortex-core elements. These two kinds of ex- 
citations are not independent but instead are hybridized 
in a way-reminiscent of the case of excitons in dielectric 
crystals Q. This analogy will guide us in order to build 
a phcnomcnological Hamiltonian. This phcnomenologi- 
cal model will allow us to derive analytical expressions 
for measurable quantities and to provide an interpreta- 
tion for the high-momentum roton region of the spectrum 
where the variational approache becomes more involved. 

The physical motivation leading to the assumption 
of the existence of localized excitations stems from the 
following remarks. In any attempt to reproduce quart, 
titatively the spectrum using a variational approachEl 
more localized structure needs to be introduced into the 
excited-state wavefunction. In addition it is generally 
accepted that the superfluid supports quantized vortices 
with microscopic coresD. In order to describe them, it is 
relevant to consider □ exchange effects associated to lo- 
calized and closed exchange rings of atoms. Therefore, 
we view the system as being a liquid at large distances 



while locally, on the scale of some interatomic distances, 
we consider localized quantum clusters playing the role 
of the vortex-core elements. Since the helium atoms are 
indistinguishable, we cannot split them into two classes 
(i.e. either embedded in a long range density mode or in a 
localized cluster) but we suppose instead that each atom 
retains both characteristics. As a consequence, these two 
kinds of excitations are coupled and the spectrum results 
from this hybridization. From a methodological-point 
of view, our approach is close to that of Hopfield □ who 
considered excitons in dielectric crystals as localized exci- 
tations whose interactions are mediated by the photons. 
Here, along the same line, the interaction between the lo- 
calized excitations is mediated by the Feynman density 
fluctuations. 

To make those qualitative remarks more precise, we 
first consider the Hamiltonian Hq = J^ fc e(/c)a/Jafc de- 
scribing the pure density modes obtained using the 
Feynman ansatz a with the energy e(fc) = 2 mS(k) > 
where m is the mass of the helium atoms and S(k) 
the structure factor of the liquid. Then, we assume 
the localized excitations to be two-level systems of en- 
ergy hujQ, such that the corresponding Hamiltonian is 

H°ioc — Ylk^ (t>k?bk + |^ . The operators a and b 
obey bosonic commutation relations and commute be- 
tween themselves, but for the b operators the bosonic 
character is only approximate and holds in the limit of 
a low density of localized modes 0. To write down the 
part describing interactions between the two sets of exci- 
fions, we follow the approach of Hopfield and Anderson 
and assume that the role of the density fluctuations 
is to induce an effective interaction between the localized 
modes which within the dipolar approximation allows to 
write the effective Hamiltonian 

H loc = Y,(huj + X{k)) (bjb k + i 

+J2 X ( k ) (ftfc^lfc + hb-k) (1) 

k 

where X(k) is a (real and negative) matrix element which 
depends on the microscopic details of the dipolar inter- 
action. The dipolar approximation is the simplest limit 
for which the resulting Hamiltonian can be analytically 
worked out. Moreover, it corresponds to the physical pic- 
ture of the localized- modes as vortex-core elements. Such 
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microscopic line-elements may be treated as local dipo- 
lar defects with respect to the surrounding fluid. As we 
shall see later, in this limit, the knowledge of the exact 
expression of X(k) is not necessary. The coupling be- 
tween the phonons and the localized modes is described 
by the Hamiltonian U 

H c = ^ (K k > u a )b k + uo)ak) {a^ + «-fc) + h.c. (2) 



where the two functions A and /x are given by A(fc) 



i?kj o(-3^)^ and fj,(k) 



f- 3X(k) 



Hi„r is di- 



agonalized using the Bogoliubov transformation (3 k — 
u{k)b k + v(k)b' -k, which the resulting spectrum 



E = ^huj {huj + 2X{k)) 



(3) 



The two functions u(k) and v(k) are given by u 2 (k) — 
+ 1} and w2(fc) = _ i). The 

Hamiltonian Ho + H°i oc + H c is quadratic and is di- 
agonalized by means of the canonical transformation 
a k = Aa k + Bb k + Crf- k +Dtf- k and a k = Ba k + Ab k + 
Da^ - k + C6^_fc. The corresponding dispersion relation 
is 



e 2 W 
E 2 {k) 
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X{k) 



hu 1 _ (Hfi)2 



(4) 



We first notice that taking the coupling X(k) between 
the two sets of modes to zero, we obtain, as expected, 
the two solutions E(k) — e(k) describing a pure den- 
sity mode and E = huiQ for the localized two-levels. A 
non zero coupling hybridizes these two sets of excita- 
tions. We take for the energy hui the highest value of the 
phonon-roton spectrum i.e. where it terminates, namely 
Hluq = 2A where the energy A corresponds to the ro- 
ton minimum. New experimental results at higher values 
of k confirm this behaviour at the termination pointEJ. 
To solve for the dispersion relation and to obtain the 
phonon part of the spectrum, we use (3) into (4). This 
gives E(k) — |e(fc) i.e. an expression independent of 
the matrix element X(k). Using now this latter expres- 
sion into (4), we obtain the other branch describing the 
hybridized local modes at E = 2hujQ = 4A which is as 
well independent of X(k). We emphasize again that this 
is a consequence of the particular choice of the dipolar 
approximationO which is manifested in the expressions 
for the coupling functions A and /i. As a result of the 
hybridization, the energy spectrum of the density fluc- 
tuations is shifted by a factor two towards lower ener- 
gies relatively to the Feynman ansatz, and the localized 
modes still have a constant energy although twice the 
original value. 

To compare our results with the experimental data 
we shall consider first two independent sets of results 
namely the measurements of the energy spectrum and of 
the structure factor S(k). We obtain from (4) the relation 
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E(k) 



fifk 



AmS(k) 



(5) 



for the lower branch which we compare to the experi- 
mental r esult s obtained-at-two different pressures both 
for £(fc)By and S{kj$M (Figs.l and 2). The main 
discrepancy is obtained in the low momentum region 
(k < lA ), i.e. below the maxon momentum. As k — > 
the experimental data is close to the Feynman result i.e. 
twice the value we obtain in (5). This signals a fail- 
ure of the dipolar approximation for long wavelengths. 
Our description becomes meaningless when the phonon- 
roton branch terminates since the Feynman spectrum ap- 
proaches that of a free particle. Around the roton mini- 
mum, and over a large range of momentum the result (5) 
provides a good fit to the experimental results. Feynman 
himself in his original paperEJ noticed the factor 2 dis- 
crepancy between his ansatz and the experimental results 
around the roton minimum. Here we deduce this factor 
from an exact analytical solution. This interpretation 
of the roton excitation as resulting from the dipolar hy- 
bridization of two separate excitations is-.to,be compared 
with the approach of Clyde and Griffin E3>U which uses a 
dielectric formalism to describe hybridized phonons and 
free-particle excitations. 

The neutron scattering intensity is a direct measure of 
the density fluctuations in the liquid and is obtained from 
the dynamical structure factor S(k,uj). It ia_.acceptcd 
since the work of Miller, Pines and Nozieres Ii3 that we 
can split S(k,cu) into two parts, S(k,uj) = N Z{k)8{hu> — 
e k ) + S , ^ 1 ^(fc,w) where the first term accounts for single 
quasi-particle excitations while the second describes mul- 
tiparticle excitations. This separation is quite easy to jus- 
tify at low temperature and low momentum, typically for 
k < 0.5 A where we obtain Z(k) = S(k), which results 
also from the Feynpian theory. The comparison with the 
experimental data EH shows that although it works in the 
low momentum regime mentionned above, it fails to de- 
scribe the behavior at higher momentum, except perhaps 
for the position of the maximum. In contrast to this, we 
do not consider here a precise decomposition of the struc- 
ture factor. Since the density fluctuations are hybridized 
with the localized modes, we expect the differential cross- 
section Z{k) for the excitation of a single quasi-particle 
to be proportional to (&fc^Lfe) = u k v k , where the ex- 
pectation value is calculated in the ground state of the 
Bogoliubov pairs. Then, Z(k) = 4irk 2 Iou k v k where Iq is 
a normalization constant. Using our previous expression 
for u k and v k we obtain u k v k — j ^gpr which together 
with (3) and (5) gives 
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(6) 



From the two independent measures of E(k) and S(k) we 
obtain using (6) a theoretical expression of the differen- 
tial cross-section which fits well the experimental results 
as shown in Fig. 3. In particlular from (6) it follows that 
the intensity of the single quasi-particle branch vanishes 
when E(k) — > 2 A which is a new result. Moreover, in 
the low momentum limit, we recover the expected pro- 
portionality between S(k) and Z{k). 
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We found a second branch of excitations at the con- 
stant energy E = 4A. It describes localized excitations 
of energies twice the bare vortex core energy. It is sugges- 
tive to interpret this mode as a single vortex loop whose 
radius can be calculated using a Feynman-type formula 
r the energy of the circulating current of a vortex-loop 
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vortex— loop 
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where at T = 0, the density of the superfluid is p = p s 
and a is the core size equal to the atomic radiusD namely 
a ~ 1AA. The radius R is obtained taking E vortex = 
4A = MAK. This gives R ~ 5.0A which is the ex- 
pected size for the smallest vortex-loop. A further ex- 
perimental evidence in favour of this interpretation is 
provided by critical velocity experiments LL3. In phase- 
slippage studies of the critical velocity through an ori- 
fice, the critical velocity is driven by the thermal nucle- 
ation of vortex-loops. The corresponding energy E v is 
determined by the nucleation rate T which, using the 
Arrhenius law V = Toexp ^j^p^j , is found El to be 

E v ~ 33 ± 5K, which is indeed very close to our re- 
sult 4A = 3AAK. Finally, the upper critical velocity v c 
may be estimated as the velocity of the vortex-loop itself 
□ namely v c = y^r}'* 1 (tt) — 20mA&, a value close to the 
largest measured critical velocity iB. 

Another relevant set of experiments we consider is pro- 
vided by the Raman scattering around k = 0. It has been 
found that besides a peak at E^j 2A, there is an addi- 
tional contribution at E = 4A Ej which we associate to 
the vortex-loop. The peak at 2A is usually interpreted as 
a two-roton excitation and therefore the additional con- 
tribution is viewed as a four-roton excitation. This inter- 
pretation suffers nevertheless from the fact that there is 
no three-roton peak. In our model, the lowest excitation 
energy of a vortex-core is given by 2A and not by A so 
that we do not expect any contribution at 3 A. The peak 
at 2A which does not appeat,in the hybridized spectrum 
will be discussed elsewhere cA. 

In conclusion, we have shown that the roton region 
of the excitation spectrum of superfluid 4 -He can be de- 
scribed quantitatively by assuming the existence of two 
kinds of excitations. One is provided by the delocal- 
ized density fluctuations given by the Feynman varia- 
tional ansatz. In contrast, the second set of excitations 
describes the short range order in the system around a 
microscopic vortex-core element. We have assumed that 
these excitations are coupled in a way reminiscent from 
the case of excitons in dielectric systems. By writing a 
phenomenological Hamiltonian together with the dipo- 
lar approximation for the coupling, we obtained an ex- 
citation spectrum which involves two sets of hybridized 
modes. One corresponds to phonons shifted by a factor 
of two towards lower energies relatively to the Feynman 
result. The second set is interpreted as the intrinsic quan- 
tized vortex-loop modes of the superfluid. This picture 
provides analytical expressions which describe quantita- 
tively a broad range of experimental results. In addi- 



tion, it may help understanding the persistence of the 
phonon peak above T\ which appears to be shifted up- 
wards in energy together with a sharp drop of both the 
maxon and roton peaks. Since the localized excitations 
depend on the superfluid order, they will vanish at the 
transition to the normal state, while the Feynman density 
modes will remain unaffected and only shifted upwards 
in energy in the abscence of dipolar coupling X{k). Our 
model provides a new picture of the nature of the roton 
excitation and may allow for a unified treatment of the 
phonon-roton elementary excitations and of the quan- 
tized vortex-loops excitations which are both unique to 
the superfluid phaseE3. 
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FIG. 2. Comparison between the experimental structure 
factor S(k) (points) and the expression (5) (solid line) for 
the same two pressures as in fig.(l) where the energy E(k) is 
obtained from independent measurements [10,11]. (a) circles 
[12], squares [14], (b) circles [13]. 
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FIG. 1. Comparison between the experimental energy 
spectrum [10,11] (points) and the theoretical expression 
(5) (solid line) where the structure factor S(k) is obtained 
from independent measurements [12,13]. (a) and (b) corre- 
spond respectively to the saturation vapor pressure and to 
P=24 atm. The dashed line at 4A indicates the position of 
the branch of the vortex-loop excitations. 
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FIG. 3. Comparison between the experimental scatter- 
ing cross-section [10] Z(k) of single quasi-particle excitations 
(points) at 1.1 K and the theoretical expression (6). The two 
curves are obtained using respectively in (6) the experimental 
results for S(k) [12] (dashed line) and for the energy E(k) [11] 
(solid line). 
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